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We show that spin-polarized current flowing through a ferromagnet leads to predominantly in- 
coherent magnetic excitations. We describe these excitations by an effective magnetic temperature 
rather than a coherent precession as in the popular spin-torque model. Our model reproduces all 
the essential features of the experiments, and gives several predictions that distinguish it from the 
spin-torque model. 



I. INTRODUCTION 

Electrons injected into a ferromagnet lead to excitation 
or reversal of magnetization. This effect was predicted by 
Slonczwesk;^ and Berger,^ and later observed experimen- 
tally.3-4 j^jQg^ ^Yic magnetization switching experiments 
are performed with trilayers of structure F1/N/F2, where 
a thick and usually extended ferromagnet Fi plays a role 
of electron current polarizer not affected by the current, 
N is a nonmagnetic spacer separating the magnetic lay- 
ers, and F2 is nanopatterned into a typically elongated 
nanopillar of submicrometer dimensions. The current- 
driven magnetization switching of the nanopillar F2 usu- 
ally occurs between two well-defined orientations along 
the magnetic easy axis defined by the nanopillar shape 
anisotropy. This is the specific experimental situation we 
will be dealing with in this paper. 

Both of the original modelsi*^ in different ways rely 
on transfer of spin from the polarized electron current 
to the ferromagnet. Berger considered generation of 
magnons by spin-flipping of electrons, driven by the spin 
accumulation. Slonczewski considered the electron spin 
component transverse to the magnetization, which is ab- 
sorbed by the ferromagnet due to a combination of spin- 
dependent reflection at the interfaces and averaging of 
the electron spin precession phases in the ferromagnet. 
The resulting torque drives the magnetic dynamics. In 
another model, proposed later by Heide et al.,^ the effect 
of current was described as a nonequilibrium exchange 
interaction between Fi and F2. This interaction was sim- 
ilar to an effective current-dependent field. Later exper- 
imental results were found to be inconsistent with this 
effective field modeliS*2i£ 

The spin-torque model (ST) has become dominant in 
the interpretation of experimentsi'^i^i'^i^i^'^i^^i^^i^^i^^i^°i^^ 
A few exceptions include some of the point con- 
tact measurements, interpreted in terms of resonant 
electron-magnon scattering, and thermally acti- 
vated switching, interpreted in terms of incoher- 
ent current-driven excitations, described as an effec- 
tive magnetic temperatureii2*ii*i24i^ With various ex- 
tensions and modifications, incorporating incomplete 
current polarization, band structure, spin accumula- 
tion, thermal activation, and inhomogeneous magneti- 
zation states, ST is also dominant in theoretical re- 
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Despite the overall success of ST, several fundamen- 
tal issues remain unresolved. First, approximation of 
the magnetization by a classical macrospin is essential 
for the coherent dynamics predicted by the model. On 
the other hand, finite-wavelength magnetic excitations 
(magnons) are also found to be excited in certain ex- 
tensions of ST, thus violating the assumptions underly- 
ing the modeliSi A large number of the excited magnetic 
modes then need to be considered to adequately describe 
the current-driven magnetic dynamics. This situation is 
fundamentally different from the magnetic excitation by 
ac field in transverse ferromagnetic resonance (FMR) ex- 
periments, where only uniform precession is excited by 
the ac magnetic field, because nonuniform excitations do 
not couple to this field. On a more fundamental level, 
the relation of ST to the quantum-mechanical picture 
involving electron-magnon scattering is not understood. 
Among related issues, energy transfer from the current 
to the magnetic system, and the magnetic dynamics in- 
duced by scattering of a single electron, are not addressed 
by the classical ST. 

The first goal of this paper is to critically examine the 
fundamentals of ST, and its relation to the quantum- 
mechanical electron-magnon scattering picture. In Sec- 
tion IHI we consider a toy- model quantum analog of ST, 
i.e. a model involving the quantum dynamics of both the 
nanopillar magnetization, treated as a macrospin, and 
current-carrying electrons. In addition to the ST result 
for the coherent magnetic excitations, this model gives in- 
coherent (i.e. not described as a classical precession) ex- 
citations. In the typical experiments, where the polariza- 
tion of the current is initially e.g. antiparallel to the mag- 
netization, the coherent excitation vanishes, while the in- 
coherent excitation rate is largest. Current-driven gener- 
ation of finite- wavelength magnons, and magnon interac- 
tions further decrease the coherence of the excitation. In 
Section UTTl we use a ballistic scattering approach, where 
the magnetic excitations are approximately described as 
increase of a magnetic temperature, to obtain the magne- 
tization switching in nanopillars. Our model reproduces 
all of the experimental current-switching behaviors, and 
gives several predictions that distinguish it from ST. 
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II. TWO-SPIN SCATTERING MODEL 

In this section, we solve a simple model for the interac- 
tion of a conduction electron with a ferromagnet F2, and 
show how the result relates both to ST— and electron- 
magnon scattering framework^ We model F2 by a large 
spin in the state |S >= |L, M >= \M >, where L is 
fixed, and M is its projection on the quantization axis 
z, aligned with the effective magnetic field H^jf- H^ff 
includes both the applied field and the anisotropy field 
of the nanopillar. We define the exchange stiffness of a 
ferromagnet as the exchange contribution to the disper- 
sion of the spin- waves in the bulk of the ferromagnet. In 
typical experiments, the ratio of the exchange stiffness of 
F2 to its largest dimension is significantly smaller than 
the typical energies of conduction electrons scattered by 
this layer (measured from the Fermi level). Thus, finite- 
wavelength magnetic modes are expected to be excited 
by the current, which cannot be described in the fixed- 
L (macrospin) approximation. When we relate our re- 
sults to experiments in Section IIIII the macrospin ap- 
proximation is not used. In the present Section, we use 
the macrospin approximation both for comparison to ST, 
which makes a similar approximation, and due to the 
transparency of the results. 

We assume that F2 is a transition metal. The magnetic 
properties of transition metals are dominated by the 3d 
electrons, while transport occurs through the more mo- 
bile states of predominantly 4s type. The hybridization 
between the conduction electrons and the magnetization 
is approximately described by a Stoner exchange poten- 
tial Vex = iJ/L)S ■ s, where s is the conduction elec- 
tron spin. The pre-factor {J/ L) explicitly accounts for 
the independence of the exchange potential of the size of 
the ferromagnet. If the orbital part of the electron wave 
function is included, Vex gives spin-dependent reflection 
at the magnetic interfaces. We omit the orbital part, 
concentrating on the spin dynamics. 

Hgf f can be ignored when solving the problem of scat- 
tering of the conduction electrons by F2 , since the associ- 
ated Zeeman energy is typically five orders of magnitude 
smaller than Vex ■ When an electron is far from the ferro- 
magnet before scattering, it is in the state |s >= a\ t> 
-|-/3| |>. To solve the scattering problem, we express the 
electron wavefunction through the angular momentum 
eigenstates 

V2LTlti = V2LTl{a\ t> i>)\L,M >= 
aVL + M + \\L + l/2\M + 1/2 > + 
/3VL-A/+l|i + l/2;M-l/2> - (1) 
a^/L - M|L - 1/2; M + 1/2 > + 
P^JL + M\L - 1/2; M- 1/2 > . 

Here \j]jz > denote the states with total spin j and z- 
projection ■ Vex = J foi' the first two terms on the right, 
and Vex — ~J^^jr ~ ^^^^ terms. We turn 

on the interaction for the time an electron at the Fermi 



energy spends in F2, t « 10~ sec for a several nm thick 
magnetic layer. In transition metal ferromagnets. Vex ~ 
— 1 eV gives a similar time scale for the spin dynamics, 
Vex/^ ~ 1 — 10 X 10~^^ sec. Scattering gives a phase shift 
« 10° - 10^ between the j = L - 1/2 and j = L+ 1/2 
terms in Eq. ^ The final state in the basis of individual 
spins is 

(2L-f 1)^^/ = 

a\ t> \M >[L + M + 1 + e"^(L - M)] + 

a\ l> \M+1> v/(i-M)(L + M + l)(l - e**) + (2) 

p\ l> \M >[L-M + 1 + e'^{L + M)] + 

p\ T> \M - 1 > v/(i + M)(^-A/ + l)(l - e'^). 

To compare with ST, we first consider L = AI, an 
example of a coherent state of F2. We calculate the elec- 
tron spin components. Since Vex conserves the total spin, 
identical results are obtained with similar calculations for 
S. Take a, (3 real, then only x- and z-components of the 
electron spin before scattering are finite 

< i^.\Sx\'^r a/3, < i^^\sM^ >= {o? - f (3) 

For the final state (Eq.EJ, we first find expectation values 
of spin components, and then average over due to its 
strong variation with t for different electron paths in F2. 
Analysis involving an actual integration over the different 
paths is expected to yield similar resultsi ^^i^° 

< V'/ISylV'/ >ai,= 0, < >av= ajS / {2L + 1) 

O T 

< V/k.lV'/ >a.^ [o? - /32[1 - ^^^^])/2. (4) 

To avoid confusion, we note that ST considers the ref- 
erence frame set by the orientation of S. This frame 
coincides with our x, y, z frame only in the special case 
M = L. The X- and z-components of electron spin are 
then usually referred to as longitudinal and transverse 
spin components, respectively, as defined with respect 
to the orientation of S. For large L, Eq. ^ shows that 
all the transverse electron spin is absorbed by S, while 
the change of the longitudinal spin is small, As^ cx 1/ L. 
If Heff = 0, S subsequently relaxes to a classical state 
with a new direction determined by the spin conservation 
sin(6') w As^/L, or « Here 9 is the angle between 
the new orientation of S and the z-axis. A similar con- 
clusion is reached by STi 

A different result is expected at finite Heff along the 
z-axis. For the sake of argument, assume that scattering 
of electrons classically tilts S away from its equilibrium 
orientation along Hef / . S subsequently starts precessing 
around Heff, relaxing back to its equilibrium state. It 
is well established in the ferromagnetic resonance (FMR) 
experiments, that the magnetic relaxation can be sepa- 
rated into two channels, characterized by the relaxation 
rates Ti and 

Ti characterizes the relaxation of the total magnetic 
energy, mostly due to scattering with the conduction elec- 
trons and phonons. In current-driven switching experi- 
ments, the magnetic energy relaxation into conduction 
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electrons is partly blocked by the same mechanism that 
gives rise to the current-driven excitation. Thus, Ti is 
expected to be larger than in the FMR experiments with 
films of the same thickness as F2 . The magnetic energy 
relaxation is measured as a variation oi n = L — M , 
which we call the number of magnons. This definition 
needs to be expanded when including nonuniform mag- 
netic excitations, which reduce the length of S. We 
do not perform the Holstein-Primakoff transformation, 
so the magnons only approximately obey Bose-Einstein 
statistics at n << L.'^^ 

T2 characterizes the relaxation of the uniform preces- 
sion, which, in addition to the processes contributing to 
Ti , also includes scattering of the uniform precession into 
other magnetic modes, through nonlinear magnetic in- 
teractions. Generally, T2 < In current-switching 
experiments, the nonlinear magnetic interactions are en- 
hanced (compared to the values for small amplitude FMR 
in bulk ferromagnets) by high excitation levels associated 
with the switchingjSl and by scattering at the surfaces of 
the thin layer ¥2^^ 

The relation T2 < Ti allows for nonequilibrium popula- 
tions of all the magnetic modes, even if only the uniform 
precession is excited by the current. Moreover, we argue 
below that the current-driven excitations themselves (re- 
gardless of the T2 processes) are at least partly incoher- 
ent, and many finite- wavelength modes are directly ex- 
cited by the current. Thus, we expect that at sufficiently 
long excitation times (e.g. when a dc current is applied), 
the uniform precession is only a small part of the gen- 
erally incoherent steady magnetization state. General 
arguments of statistical physics tell us that the nonlinear 
magnetic interactions tend to thermalize the magnetic 
excitations. In Section IIIII the excited state is approx- 
imately described by an effective magnetic temperature 
T^(/)i2iiiiiM2i& Tm{I) characterizes the energy stored 
in the magnetic system, and gives approximate popula- 
tions of the magnons n,- « — ,„ ,} 1 — 7. Here Ei are 
magnon energies. 

We estimate the total rate of magnon generation in the 
macrospin approximation for F2 . Eq. 0] gives the trans- 
verse spin absorbed by S — — As^ ~ —j-. If one 
assumes that that the longitudinal spin-transfer is asso- 
ciated with coherent magnon generation, the correspond- 
ing initial angle 6 of classical precession of S around Hef / 
would follow from cos(6') w 1 — Asz/L, giving 6* « ^. 
It is always larger than the estimate 6* « ^ based on 
the transfer of the transverse spin (see above). This 
means that scattering from a coherent state of S (e.g. 
\L,L >) to another (tilted) coherent state constitutes 
only a fraction of all the scattering processes. The trans- 
fer of longitudinal spin is largest when a = 0; then the 
excitation is completely incoherent, since Asx = 0. Inci- 
dentally, in typical current-driven switching experiments, 
the magnetizations Mi and M2 of layers Fi and F2, cor- 
respondingly, are either parallel or antiparallel to each 
other liSiiiiSiiSiii M2 is related to S by the gyromagnetic 



ratio. Electrons flowing through F2 then do not have a 
spin component perpendicular to S. In this case a/3 = 0, 
and the transverse spin transfer, and thus coherent exci- 
tation, vanishes. 

Eq. 2] clarifies a popular misconception that the trans- 
verse electron spin component can be elastically (without 
electron spin-flipping) transferred to the magnetization. 
In this context, the probability of electron spin-fiipping 
is identified with the longitudinal spin transfer. In the 
framework of the macrospin model, at large L the prob- 
ability of electron spin-flip (or equivalently magnon gen- 
eration) decreases as 1/L. On the other hand, the x,y 
(transverse at small excitation levels) components of an- 
gular momentum, associated with excitation of uniform 
precession, are proportional to L (per magnon). Thus, in 
the large-L limit, the coherent magnon generation gives 
transverse spin tranfer (per scattering electron) indepen- 
dent of L, consistent both with ST and the formalism of 
electron-magnon scattering. These arguments, of course, 
are valid only in the macrospin approximation (see be- 
low). 

Because we do not consider the spatial components of 
the magnetic system and the electron wave function, our 
model does not apply to finite- wavelength magnons. We 
can qualitatively describe such excitations by breaking up 
the magnetic layer into N sufficiently small interacting 
elements with spins Li — L/N,i — 1..N, in the spirit of 
micromagnetic models. The probability to excite one of 
these elements and not the others is then a rough measure 
of the probability to excite a magnon with wavelength 
of about the size of the element. According to Eq. 01 
Asz oc so the probability to excite magnons with 

shorter wavelengths may be expected to be higher. 

In the limit of small exchange stiffness of F2, we can 
replace our macrospin S pointing along the z-direction, 
with a system of 2L weakly interacting electrons with 
spin Sz = 1/2. Our problem of conduction electron scat- 
tering on a macrospin can now be replaced with another, 
where we add one electron with a spin s, and then remove 
one electron from the system. When the conduction elec- 
tron spends a sufficiently long time in F2, Vex > fi/t, the 
outgoing electron will be any one of the 2L+1 electrons in 
the system, and will have spin projection Sz = 1/2 with a 
large probability P > 1 — 1/ (2L). The difference between 
this result and Eq. 2|is due to the excitation of magnons 
with finite wavelengths. This (grossly oversimplified) ar- 
gument demonstrates again that all the magnetic modes 
are excited by the current, limited only by energy con- 
servation. It can be estimated that in typical current- 
switching experiments at not too large applied fields, en- 
ergy conservation can be satisfied for a large number of 
magnetic modes: typical energies and spacings between 
low-energy magnetic modes in transition metal magnetic 
nanopillars of submicron size are ~ 10/ieV, while typical 
thermal electron energies are « 25 meV at 295 K, the 
temperature at which most experiments are performed, 
or w 0.3 meV at 4.2 K. Moreover, at typical experimen- 
tal switching currents of several mA, the electron energies 
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due to the voltage applied to the sample are « 1 meV, 
also much larger than the energies of the lowest magnetic 
modes. 

Relating our arguments for current-driven excitation 
of finite- wavelength magnons to ST, we notice that the 
finite-wavelength magnons do not conserve the length of 
S, thus violating the coherent approximation underlying 
ST: generally, Asx — for these modes: when they are 
excited by the current, spin transfer transverse to the 
total magnetic moment of F2 does not occur. Recently, 
local arguments in the framework of ST were made in fa- 
vor of excitation of finite-wavelength magnetic modesn^i 
Just like the original ST dealing only with the uniform 
precession, this model captures only the classical (i.e. in 
this case locally coherent) magnetic dynamics. The local 
spin-torque effect also requires a semiclassical approxima- 
tion for the conduction electrons, in contrast to the orig- 
inal ST, which treats electrons quantum- mechanically.^ 

Wc now consider a more general case of scattering with 
S initially in \M >,M < L state. From Eq.[21 
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2(2i + l)2 



[|/3|2((2i-Hl)2-(2Af-l)2)- 
|a|2((2L + l)2-(2Af + 1)2)]. (5) 



The first term on the right describes the probability of 
magnon generation by a spin-down electron, the second 
describes magnon absorption by a spin-up electron. An 
important limit of Eq. is n << L,M. The magnon 
emission probability is then 



P. 



\{2L + If - {2M - If 
2(2L+ 1)2 



,1 



L 



(6) 



giving both spontaneous (independent of n) and stimu- 
lated (cx n) generation of magnons by the currenti^ The 
magnon absorption probability is 



Pn 



a 



^((2i-|-l)2-(2Af + 1) 



2{2L + 1Y 



L' 



(7) 



Eqs. I6I7I reproduce the general Einstein formulae for 
bosons excited by an external field, valid also for the 
generation of finite-wavelength magnons. Although the 
Einstein formulae give correct values for magnon gener- 
ation rates, they do not contain information about the 
coherence of these excitations. ST shows that coherence 
is related to the transverse spin-conservation, not ade- 
quately addressed in the formalism of electron-magnon 
scattering. 

Another important limit of Eq. (Sjis < M « L, at 
high excitation level of S. This limit likely corresponds 
to the high-field reversible switching regime in the exper- 
iments with nanopillarsi^. For < M << L, Eq. [3 gives 



i«n/2, 



(8) 



expressing saturation of the magnon excitation probabil- 
ity. The z-component of the electron spin, which is now 
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FIG. 1: Schematic of our model for the spin-dependent trans- 
port through the nanopillars, as explained in the text. 



perpendicular S, is completely absorbed by S, in agree- 
ment with ST. However, even in this limit, the excited 
state of S can be completely incoherent: all the projec- 
tions (x, y, and z) of S in the state |L, > vanish, i.e. it 
is a purely quantum state. Therefore, Eq. |H1 represents a 
generalization of the ST result to the highly excited in- 
coherent magnetization states. Below, we argue that the 
saturation of magnon generation is likely responsible for 
the behaviors associated with the reversible switching. 

To summarize this section, Eq. 0] shows that ST cap- 
tures only a small, coherent part of the current-driven 
excitations. The Einstein formulae Eas. 16171 describe the 
current-driven excitations in terms of electron-magnon 
scattering, but they do not provide information about 
the level of coherence of these excitations. We gave ar- 
guments that coherence may not be important in typ- 
ical experiments. With this assumption, in Section IIIII 
we use the electron-magnon scattering formalism, and 
effective temperature description for the magnetic ex- 
citations, to reproduce the experimental current-driven 
magnetization switching behaviors. 



III. 



CURRENT-DRIVEN SWITCHING 



In this Section, we show how incoherent current-driven 
magnetic excitations of F2 , rather than a coherent rota- 
tion of S due to the spin-torque mechanism, may lead 
to the experimental current-driven behaviors. To com- 
pletely eliminate the coherence associated with the spin- 
torque effect, we assume that equilibrium orientation of 
the magnetization Mi is (anti)parallel to M2, due to the 
shape anisotropy of the nanopillar F2. 

Evolution of a system with many excited modes can 
be described in terms of the density matrix dynamics. 
Besides analytical and computational difficulties, such a 
description would involve a large number of presently un- 
known parameters. Instead, we now show how an approx- 
imate effective temperature description of the magnetic 
system excited by current, and current-driven magnon 
generation rates given by Eqs. I6I7I reproduces typi- 
cal current-driven switching experiments with nanopil- 
^g^j.gj^ji3ji5j^ji7 rpj^jg analysis does not prove that the 
current-driven excitations are thermalized. Rather, it 
shows that coherence of the excitations is not essential 
in these experiments. 
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Fig. ^ illustrates our quasi-ballistic model of transport 
through a F1/N/F2 trilayer, similar to the circuit model 
of Brataas et al^ The nanopillar F2 is represented by 
a high-resistivity constriction, separating low-resistivity 
reservoirs. If we neglect inelastic scattering in F2 (includ- 
ing spin-flipping), and separately consider two different 
spin channels, this system can be treated in the Landauer 
formalism. To account for two separate spin channels, 
one can formally replace each reservoir in Fig. □ with 
two completely spin-polarized reservoirs. The electron 
spin-flipping can then be treated as transmission between 
spin-up and spin-down reservoirs. 

We model the spin-polarizing effect of the large ferro- 
magnet Fi by different spin-up and spin down potentials 
in the left reservoir, V'^ and VK The spin states are 
defined with respect to the majority spin orientation in 
F2 (i.e. opposite of the magnetization M2). We take 

V = in the right reservoir, i.e. neglect the spin accu- 
mulation there. l\ are spin-up and spin-down cur- 
rents through the constriction. Positive / is from Fi to 
F2. We define AV, p by AV ^ - ^ 2pV. Here 

V !v (V'^ + V^)/2 is related to the voltage across the 
sample. If F2 is removed, p becomes the current polar- 
ization in the constriction. Spin-dependent scattering at 
each of the interfaces of F2 is described by conductances 

= 2iyV'^, = We introduce an addi- 

tional ballast resistance i?o in both spin channels, giving 
a correction to these expressions. Rq accounts for contri- 
butions from Fi, and other bulk and contact resistances. 
It controls the asymmetry of the switching currents in 
the opposite directions. We neglect multiple scattering, 
so the interface resistances add. 

At small H, the current-driven switching occurs 
through thermal activation at relatively low magnetic 
excitation levels, with magnon generation rates approx- 
imately described by Eqs. 16171 ^'^^ A diffusive model for 
the stimulated current-driven magnon emission, driven 
by spin-accumulation, was developed in Ref. 2, where 
the energy was used as the only parameter character- 
izing the electron distributions. The ballistic model of 
Fig. ^ needs to incorporate the nonequilibrium distribu- 
tions in momentum space, and hot electron transport 
through the constriction. However, such extensions of 
Ref. 2 would have limited physical meaning, since the 
semiclassical approximation, involved in the distribution 
function formalism, breaks down on the relevant length 
scales of the constriction. 

Instead, we use ballistic scattering theory arguments 
to obtain an expression for the spin-flip scattering rates 
in the constriction. For simplicity we neglect magnon 
energies, which (for the low-energy modes) are smaller 
than the typical conduction electron energies at the bias 
voltages leading to switching, and consider spin-flipping 
of electrons both transmitted and reflected at the inter- 
faces. To avoid complications associated with scatter- 
ing among different spin-channels, we separately consider 
spontaneous and stimulated magnon emission. 

First, we assume p = 0, i.e. the current from Fi is un- 



polarized. In this limiting case, the spin-flip scattering 
back into the same reservoir is suppressed. Both spin- 
up and spin-down electrons experience spin-flip scatter- 
ing with identical matrix elements, as follows from the 
unitarity of the scattering matrix. Because of the spon- 
taneous emission (see Eq. O , the total rate of magnon 
generation by the spin-down electrons is higher than the 
magnon absorption rate by the spin-up electrons, and 
is proportional to the number of electrons transmitted 
by F2, Tsp « A|ey|. Here A is a constant determined 
by the parameters of the constriction and the electron- 
magnon scattering matrix element, e — — |e| is the elec- 
tron charge, and the number of transmitted electrons is 
proportional to the voltage drop across the constriction 
(Ohms law). The spin-polarization due to F2 is usually 
small because of the ballast resistance Rq , so we do not 
consider its effect on the spontaneous magnon generation. 

We now consider the effect of finite polarization p 7^ 
in the left reservoir (Fig. essential for stimulated 
magnon generation. For simplicity, we assume here that 
every electron spin-flipping is associated with magnon 
emission or absorption. The spin-flipping due to spin- 
orbit interaction does not qualitatively change our argu- 
ments. The number of open spin-flip channels for trans- 
mission from spin-down states in the left reservoir to spin- 
up states in the right reservoir is proportional to V^, 
therefore the corresponding stimulated magnon emission 
rate in F2 is T^m ~ BieksTmV^ , with a constant Bi. 
In terms of Landauer-Buttiker formalism, BieksTm is 
related to the coefficient of transmission from the spin- 
down to the spin-up channels. We use rii « ksTm/Ei 
for degenerate magnon populations. We note that in 
the effective temperature terminology, magnon emission 
leads to magnetic heating, while magnon absorption leads 
to current-driven magnetic cooling. By a similar argu- 
ment. Tab ~ BiekBTmV^ is the stimulated magnon ab- 
sorption rate, due to spin-up electrons in the left reser- 
voir transmitted into the right reservoir with a spin- 
flip. In the remaining spin-flip scattering processes, the 
spin-up (spin-down) electrons in the left reservoir are 
reflected back into the same reservoir, with a spin-flip. 
The rate for these processes is B2ekBTm{V^ — V^). The 
total stimulated magnon generation rate is then Tst ~ 
BekBT^iV^-V^) = -2BepVkBT^, where B = B1+B2. 

In the relaxation time approximation, the spin-lattice 
relaxation rate is Tph « 7^3 (T^ — Tph), where 7 is re- 
lated to the Gilbert damping parameter in the classical 
dynamic equations.^ In thin fllms, damping is expected 
to have a large contribution from clcctron-magnon scat- 
tering.'*^ In the magnetic system with a nonequilibrium 
electron distribution, damping then becomes dependent 
on the current amplitude and polarization. We neglect 
this dependence here. In a steady state, T^t + ^sp = Tp/i, 
giving 



^ \eV\A/2_+kBTph 
^ '"^ l + 2peVB/j ■ 



(9) 



We use Eq. |51 to describe current-driven thermally acti- 
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FIG. 2: Calculated magnetization stability diagrams: (a) for 
Co/Cu/Co nanopillars at 295 K, = 2, Bj^ = 0.1 meV"\ 
(b) same as (a), but with ^ = 0, (c) for Py/Cu/Py nanopillars 
at 4.2 K, yl/7 = 0.1, B/7 = 0.3 meV^V (AP)P denotes 
a region where (anti)parallel configuration of the magnetic 
layers is stable, P/AP is a bistable region. Both orientations 
are unstable in U. 



vated switching over a magnetic barrier ?7, given by the 
shape anisotropy of the nanopillar. The switching oc- 



curs at fcsT,,,, « ;n(t^'"^oj ~ Yg, for a data acquisition 
time tex-p of 1 second per point, and an effective attempt 
rate « 10''s~"'^i^ When thermal activation is included 
in ST, the effect of / is equivalent to renormalization of 
[/f24 This can also be formally rewritten as renormaliza- 
tion of T^hr giving an expression similar to Eq. El but 
without the spontaneous term. This similarity is super- 
ficial, since ST result is based on a completely different, 
coherent picture of the current-driven magnetic dynam- 
ics. 

We use an approximation U{H,Tm) ~ 
Uo^l~T,n/Tc{l T H/Hcfm. Here the sign de- 
pends on the mutual orientations of H and M2, Tc is 
the Curie temperature of F2, and is determined by 
the shape anisotropy of the nanopillar. For a given 
value of H, we solve the equation U{H,T,n) = 16fcBT„j 
for Tjn, to obtain the effective temperature Tq{H) 
at which the switching occurs. Taking into account 
the relationships between the currents and voltages 
]/T(i) = /T(i)(i?o + 2/G'T(^)), from Eq. H we find the 
switching current 



fcB(To(F)-Tp^)(^ 



i-i-P 



-Ro- 



sign{p)eA/^ + 2epkBTo{H) B 



(10) 



For Co/Cu/Co or Py/Cu/Py (Py is typically NigoFezo) 
nanopillars with areas « we use = 30 S, 

= 6 S, p ~ ±0.6,^'* with the sign given by the mutual 
orientations of the magnetic layers. To compare with 
the pubhshed data,*^'i2-i3-i6 use ^ 1400 K and 
Uo ~ f-5 eV for for 2.5 nm thick Co nanopillars with 
typical dimensions of 140 x 70 nm, Tc = 800 K and 
Uq ~ 0.7 eV for 6 nm thick Py nanopillars with simi- 
lar dimensions. Fig.|2Ia) shows a magnetization stability 
diagram obtained from Eq. ^] for Co/Cu/Co nanopil- 



lars at T, 



ph 



295 K. A/j ^ 2, B/j = 0.1 meV" 



i?o w 1 n give a good overall agreement with the pub- 
lished data. Microscopic evaluation of A, B is outside 



the scope of this work. Within our model, A and B are 
determined by the same electron-magnon scattering ma- 
trix elements. The stimulated scattering rates are pro- 
portional to ksTmB — Brii/Ei, while the spontaneous 
emission rate is simply proportional to A with a simi- 
lar pre- factor. Thus, A/B w 20 nieV is a crude mea- 
sure of the typical energies Ei of magnons excited by 
the current. Even assuming that this is an overestimate, 
Eqs. I9I10I may be altered if the high magnon energies 
are taken into account when electron-magnon scattering 
is considered. The corresponding magnon wavelength is 
w 1 nm, supporting our estimates that the current excites 
predominantly short- wavelength magnons. 

For qualitative comparison with ST,*'^^ Fig. |2Ib) 
shows the same calculation as in Fig. |2Ia), but without 
the spontaneous term. The main difference is the lack 
of a knee in the stability line at / = 0, evident in some 
published data.^'^^ A knee at / < in some ST-based 
calculations-'^^ is different, it is due to a singularity asso- 
ciated only with the Tph = limit. The knee in Fig.[2Ia) 
is due to the contribution of the spontaneous emission in 
Eq. El A competition between the stimulated magnetic 
cooling and spontaneous heating at / > in the state 
with antiparallel magnetizations gives overall weak cool- 
ing. In some cases,"*^ the switching fields decrease both at 
/ > and / < 0, indicating that heating occurs instead 
of cooling. Spontaneous magnon emission may also be 
responsible for some of the current-driven phenomena in 
single magnetic layers As our ballistic analysis 
shows, the stimulated magnon emission vanishes in this 
case {p = 0). 

We were able to model the published 295 K Py/Cu/Py 
data with A, B values similar to those used for Co. 
Fig. |2Ic) shows a calculation for Py/Cu/Py at Tph — 
4.2 K. It reproduces the nearly "square" switching dia- 
gram. ^1 A/^ = 0.1, B/j = 0.3 meV~^, give good over- 
all agreement with the experiment. The differences with 
the 295 K values reflect a decrease of the average en- 
ergy of magnons excited by /. The dependence of A,B 
on Tph does not appear in our model, and warrants a 
more detailed analysis of the temperature dependence of 
electron-magnon scattering. 

Beyond the switching diagrams, we identify the thresh- 
old It for a large increase of Tm in the state with parallel 
magnetizations of Fi, F2 at / > 0,^^ with the diver- 
gence of Eq.Elat —eV 7/ {2pB). As we showed above, 
magnon generation saturates at large magnetic excitation 
amplitude, becoming nearly independent of Tm- Eq. (Sjis 
then replaced by a linear relation Tm ~ Tph+KpV , where 
if is a constant determined by the magnetic relaxation 
rate. A similar (with a different constant term) linear re- 
lation has been empirically proposed to explain the lack 
of temperature dependence of telegraph noise dwell times 
variations with current. ^■^ In contrast, a formal relation 
Tm oc Tph expected for the ST'^'^ does not fit those data. 
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IV. SUMMARY 

In Section^ we demonstrated that spin-polarized cur- 
rent flowing through a ferromagnet leads to both coher- 
ent and incoherent magnetic excitations (not conserving 
the expectation value of the magnetic moment). Only 
the coherent excitations are captured by the spin-torque 
model (ST). On the other hand, the electron-magnon 
scattering framework describes all the excitations, but 
does not give information about the level of excitation co- 
herence. We argue that coherence of excitations may not 
be significant or important in typical current-switching 
experiments, due to direct generation of incoherent exci- 
tations by the current, and nonlinear magnetic interac- 
tions. 

In Section ITTIl we describe an incoherent excitation 
of F2 in terms of a current-dependent magnetic tem- 
perature, and show that this approximation is consis- 
tent with the results of typical current-switching exper- 
iments. Without the spontaneous magnon generation 
term, Eq. |^ describing the effect of current on the mag- 
netization, is mathematically similar to the finite tem- 
perature result of ST. The two models, however, are 
based on completely different underlying physical pic- 
tures. Here we give two examples which distinguish 
these models. First, the spontaneous magnon generation 
term in Eq. I^is absent in ST. It gives: i) excitations at 
small Tph for layers with parallel static magnetizations, 
a nominal configuration in typical current-switching ex- 



periments; ii) current-driven excitations by unpolarized 
current )^4i2*i2ii^ iii) a knee in the magnetization stabil- 
ity curve at / = OfiSiil related to stronger current-driven 
magnetic heating than the cooling achieved by reversal 
of the current. Second, Eq. |51 suggests that the current- 
driven excitation rate may be enhanced by increasing 
the polarization^'* and electron-magnon scattering rates. 
Larger excitation levels give smaller switching currents, 
needed for possible technological applications of direct 
current-driven switching in magnetic memory devices. 
Higher scattering rates may be achieved by decreasing 
the effective magnetic volume excited by the current 
(value of L in Eq. O, through interface roughness, al- 
loying, or building complex magnetic layers with a small 
exchange strength at the interfaces. The enhancement 
of electron-magnon scattering at rough interfaces, and 
in magnetically disordered alloys, can be understood not 
only in terms of the oversimplified model of Section ^1 
but also from a more rigorous calculation of electron- 
magnon scattering amplitudes In contrast, ST gives 
reduced tranverse spin transfer for smaller magnetic vol- 
ume (oc L in Eq. 0}, and independent of the interface 
roughness and magnetic disorder. 

The author acknowledges helpful discussions with V.I. 
Kozub, N. Sinitsyn, Ya. Bazaliy, J. Bass, N.O. Birge, 
W.P. Pratt Jr., A.H. MacDonald, D.C. Ralph, A. Pert, 
support from the MSU CFMR, GSM, the NSF through 
Grants DMR 02-02476, 98-09688, and 00-98803, and Sea- 
gate Technology. 



* J. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996). 
^ L. Berger, Phys. Rev. B 54, 9353 (1996). 

^ E.B. Myers, D.C. Ralph, J.A. Katine, R.N.Louie, R.A. 

Buhrman, Science 285, 867 (1999). 
" M. Tsoi, A.G.M. Jansen, J. Bass, W.C. Chiang, M. Seek, 

V. Tsoi, and P. Wyder, Phys. Rev. Lett. 80, 4281 (1998); 

81, 493(E) (1998). 
^ C. Heide, P. E. Zilberman, and R. J. Elliott, Phys. Rev. B 

63, 064424 (2001). 
^ E. B. Myers, F.J. Albert, J.C. Sankey, E. Bonet, R.A. 

Buhrman, and D.C. Ralph, Phys. Rev. Lett. 89, 196801 

(2002). 

F. J. Albert, N.C. Emley, E.B. Myers, D.C. Ralph, and 
R.A. Buhrman, Phys. Rev. Lett. 89, 226802 (2002). 

* R.H. Koch, J.A. Katine, and J.Z. Sun, preprint. 

^ M. Tsoi, V. Tsoi, J. Bass, A.G.M. Jansen, and P. Wyder, 

, Phys. Rev. Lett. 89, 246803 (2002). 
^° J. E. Wegrowe, H. Hoffer, Ph. Guittienne, A. Fabian, L. 

Gravier, T. Wade, and J. P. Ansermet, J. Appl. Phys. 91, 

6806 (2002). 
^1 J.E. Wegrowe, lcond-mat/0306103";2003). 
^2 S. Urazhdin, H. Kurt, W.P. Pratt Jr., and J. Bass, Appl. 

Phys. Lett. 83, 114 (2003). 
" S. Urazhdin, N.O. Birge, W.P. Pratt Jr., and J. Bass, Phys. 

Rev. Lett. 91, 146803 (2003). 
" S. Urazhdin, N.O. Birge, W.P. Pratt Jr., and J. Bass, 

|cond-mat/0309191 (2003). 

J.A. Katine, F.J. Albert, R.A. Buhrman, E.B. Myers and 



D.C. Ralph, Phys. Rev. Lett. 84, 3149 (2000). 

F.J. Albert, J.A. Katine, R.A. Buhrman, and D.C. Ralph, 

Appl. Phys. Lett. 77, 3809 (2000). 

J. GroUier, V. Cros, A. Hamzic, J.M. George, H. Jaffres, 

A. Fert, G. Faini, J.B. Youssef, and H. Legall, Appl. Phys. 
Lett. 78, 3663 (2001); ibid, Phys. Rev. B 67, 174402 
(2003). 

J. Z. Sun, D.J. Monsma, T.S. Kuan, M.J. Rooks, D.W. 
Abraham, B. Oezyilmaz, A.D. Kent, and R.H. Koch, J. 
Appl. Phys. 93, 6859 (2003). 

Y. Ji, C.L. Chien, and M.D. Stiles, Phys. Rev. Lett. 90, 
106601 (2003). 

2° S.L Kiselev, J.C. Sankey, I.N. Krivorotov, N.C. Emley, R.J. 
Schoelkopf, R.A. Buhrman, D.C. Ralph, Nature 425, 380 
(2003). 

B. Oezyilmaz, A.D. Kent, D. Monsma, J.Z. Sun, M.J. 
Rooks, and R.H. Koch, Phys. Rev. Lett. 91, 067203 (2003). 
Y.B. Bazaliy, B.A. Jones, S.C. Zhang,Phys. Rev. B 57, 
R3213 (1998). 

J. Slonczewski, J. Magn. Magn. Mater. 195, L261 (1999). 
X. Waintal, E.B. Myers, PW. Brouwer, and D.C. Ralph, 
Phys. Rev. B 62, 12317 (2000). 
J. Z. Sun, Phys. Rev. B 62, 570 (2000). 
'^^ J. Slonczewski, J. Magn. Magn. Mater. 247, 324 (2002). 
M. D. Stiles and A. Zangwill, J. Appl. Phys. 247, 324 
(2002). 

M. D. Stiles and A. Zangwill, Phys. Rev. B 66, 014407 
(2002). 



8 



^® A. Kovalev, A. Brataas, and G. E. W. Bauer, Phys. Rev. 

B 66, 224424 (2002). 
^° K. Xia, P.J. Kelly, G.E.W. Bauer, A. Brataas, and I. 

Turek, Phys. Rev. B 65, 220401 (2002). 

S. Zhang, P. M. Levy, A. Pert, Phys. Rev. Lett. 88, 236601 

(2002) . 

A. Shpiro, P. M. Levy, S. Zhang, Phys. Rev. B 67, 104430 

(2003) . 

Z. Li and S. Zhang, "cond-mat /0302339" 
M.L. Polianski and P.W. Brouwer, cond-mat/0304069 
M. Sparks, Ferromagnetic- Relaxation Theory, McGraw- 
Hill, New York, 1964. 

T. Holstein and H.Primakoff, Phys. Rev. 58, 1098 (1940). 
H.J. Suhl, J. Phys. Chem. Solids 1, 209 (1957). 
R.D. McMichael, D.J. Twisselmann, and A. Kunz, Phys. 
Rev. Lett. 90, 227601 (2003). 
*® S.J.C.H. Theeuwen, J. Caro, K.P. Wellock, S. Radelaar, 
C.H. Marrows, B.J. Hickey, and V.I. Kozub, Appl. Phys. 
Lett. 75, 3677 (1999). 



J. Hohfeld, E. Matthias, R. Knorren, and K.H. Benne- 

mann, Phys. Rev. Lett. 78, 4861 (1997). 

A. Brataas, Yu.V. Nazarov, and G.E.W. Bauer, Phys. Rev. 

Lett. 84, 2481 (2000). 
''^ Y. Tserkovnyak, A. Brataas, G.E.W. Bauer, Phys. Rev. 

Lett. 88, 117601 (2002). 
'^^ R. H. Koch, G. Grinstein, G.A. Keefe, Yu Lu, P.L. Trouil- 

loud, W.J. Gallagher, and S.S.P. Parkin, Phys. Rev. Lett. 

84, 5419 (2000). 

W. P. Pratt Jr., S.D. Steenwyk, S.Y. Chiang, A.C. Scae- 
fer, R. Loloee, and J. Bass, IEEE Trans. Magn. 33, 3505 
(1997). 

S. Urazhdin, N.O. Birge, W.P. Pratt Jr., and J. Bass, un- 
published. 
*® V.I. Kozub, private communications. 

S. Urazhdin, N.O. Birge, W.P. Pratt Jr., and J. Bass, 
cond-mat/0310472 



